Solution of the multi-reggeon compound state problem in multi-co 
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We study the properties of the colour-singlet compound states of reggeized gluons in multi-colour 
QCD using their relation with noncompact two-dimensional Heisenberg spin magnets. Applying 
the methods of integrable models, we calculate their spectrum and discuss the application of the 
obtained results to high-energy asymptotics of the scattering amplitudes in perturbative QCD. 
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1. High-energy asymptotics of the scattering am- 
plitudes in perturbative QCD are dominated by the 
contribution of an infinite number of soft gluons 
exchanged in the t— channel between the scattered 
particles. Due to the remarkable property of gluon 
reggeization, it becomes possible to resum the corre- 
sponding Feynman diagrams by formulating an effective 
theory, in which collective gluonic excitations - the 
reggeized gluons, or reggeons, plays the role of new 
elementary fields 0]. In this way, the scattering ampli- 
tudes describe the propagation of the reggeons in the 
t— channel and their interaction with each other. Each 
diagram in the effective theory becomes equivalent to an 
infinite sum of the Feynman diagrams involving "bare" 
gluons. To the leading logarithmic approximation, 
a s lns = fixed as s — > oo, the dominant contribution 
comes from the diagrams with only two reggeons 
exchanged in the <— channel. It leads to rising of the 
scattering amplitudes A(s,t) with the energy s - the 
BFKL Pomeron |jj , but violates the unitarity constraints 
at high energies. To restore the unitarity of A(s,t) in 
the direct (s— , t— and u— ) channels, one has to take 
into account the t— channel exchange of an arbitrary 
number (N = 3,4,...) of reggeons. For fixed N. their 
contribution to A(s, t) satisfies the BKP equation ||. Its 
solution defines the colour-singlet compound states built 
from TV reggeons and satisfying the Schrodinger equation 

Hn^{zi,z 2 , ...,z n ) = E N ^(zi,z 2 ,...,z N ) (1) 

with the effective QCD Hamiltonian Hn acting on two- 
dimensional transverse coordinates of reggeons, Zk (k — 
1,...,N) and their colour SU(N C ) charges. In the large 
N c limit the colour structure of the Hamiltonian simpli- 
fies significantly leading to Hn — Efcli ^(zk, z*k+i) with 
Zn+i = Z\ and H(zk, Zfc+i) being the BFKL kernel Q. 

High-energy asymptotics of the scattering ampli- 
tude due to the t— channel exchange of the compound 
N— reggeon states is given by 

°° -a s E N /4 

A( S , t) ~ -i S (^s) N f - 1 UMv tiAmMt) (2) 

^ 2 (a s In a) V«»w 

with a s = a s N c /n, En being the energy of the ground 



state for (|l|) and the exponent d/v parameterizing the dis- 
tribution of the energy levels next to the ground state (see 
Eq. ( |27| ) below). The residue factors £1(2), jv(0 measure 
the overlap of ^(zi,z*2,..-,Zn) with the wave functions 
of two scattered particles. At N = 2 and N = 3 the 
corresponding compound states are known as the BFKL 
Pomeron jlj and the Odderon ]|, |i| , respectively. 

In this letter, we shall present the results of our calcula- 
tions of the ground state energy En for higher reggeized 
gluon compound states in multi-colour QCD and discuss 
their physical interpretation. 

2. Our approach to solving the Schrodinger equation 
(|lf) is based on the remarkable properties of integrability 
of the effective QCD Hamiltonian TCn in the multi- 
colour limit H [|. It turns out that, for iV c — > 00, the 
Schrodinger equation (|l|) possesses the hidden set of 
the conserved charges [Ti.N,Qk] = [<lk->Qm\ = and is 
completely integrable. It defines a quantum-mechanical 
model of N interacting particles on the two-dimensional 
plane, which can be identified as two-dimensional non- 
compact Heiseinberg spin chain model Q. The number 
of sites of the model is equal to the number of reggeons, 
N, and the spins attached to each site, S^'° and «S^'°, 
are six generators of the principal scries of the SL(2, C) 
group. Using the (anti)holomorphic coordinates on the 
two-dimensional (x, y)— plane, z = x + iy and z = x — iy, 
one represents the spin operators as the differential 
operators S° = zd z + j, S~ = -d z , S+ = z 2 d z + 2jz. 
The operators S '° are given by similar expressions 
with z, j replaced by z, j, respectively. The pair 
of complex parameters specifies the SX(2,C) 

representation. For the principal series of the 5L(2,C) 
they are parameterized by integer rij and real Vj as 
j = (1 + tij)/2 + ivj and j = 1 — j*. To match the 
reggeon Hamiltonian, Eq. (|l|), the spins have to be equal 
to j = and j = 1. In this representation, the BFKL 
kernel defining the interaction between two reggeons 
describes the interaction between the nearest spins 

JV 

H N = J2 [H(Jk.k+i) + H(J k ,k+i)] , (3) 
fc=i 

where H(J) = ip(l - J) + ip(J) - 2tp(l) with ip(x) = 
dlnT(x)/dx, Jk,k+i is the sum of two SX(2,C) spins, 
Jk,k+i(Jk,k+i — 1) = (Sk + Sfc+i) 2 with Jn,n+i = Jna, 
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and similar for Jk,k+i- Eq. ([}]) implies that the dynamics 
in the z— and z— sectors is independent on each other. 
The interaction between two sectors is ensured by the 
condition that the eigenstates VP (if, £2, %n) have to 
be single-valued functions on the two-dimensional plane. 

The two-dimensional noncompact spin magnet defined 
by Eqs. (|f|) and (|J) can be solved exactly within the 
R— matrix approach by using the method of the Baxter 
Q— operator. As was shown in the Hamiltonian is 
expressed in terms of a single operator Q(it, u) 

U N = i-^-ln{u 2N [Q(u-i(l- j),u-i(l-j)tf 

xQ(u + i(l-j),u + i(l-j))\ (4) 

and the Schrodinger equation (Q) turns out to be equiv- 
alent to the eigenproblem for the Baxter operator. Its 
eigenvalues, Q(u,u), satisfy the following conditions [Q. 
Firstly, Q{u,u) has to fulfil the chiral Baxter equation 

t N (u) Q(u, u) = (u + ij) N Q(u + i, u) + (u — ij) N Q(u — i, u) 

(5) 



withijv(u) = 2u" +q 2 u N ~ 2 + ...+q N and q = (q%, . . . , gjv) 
being the eigenvalues of the holomorphic integrals of mo- 
tion. Q{u 1 u) obeys similar equation in the u— sector with 
j = 1 — j* and q~k = The lowest integral of motion, 
q 2 , is related to the total SL(2, C) spin of the system, h, 

q 2 = -h{h-l) + Nj(j-l), h={l + n h )/2 + iv h (6) 

with rih integer and Vh real. Secondly, at u = A — in/2 
and u = A + in/2 with n integer, Q(u,u) has to be 
an analytic function on the complex A— plane except 
the points {u^ = ±z(j — m) , w^j = ±i(J — m)}, with 
m,m positive integer, where it has iVth order poles. 
The behaviour of Q(u,u) at the vicinity of the pole at 
m = m = 1 can be parameterized as ]?], |8|, ^] 

I iE ± (q) 



Q(uf 



e,uf 



C N-1 



(7) 



The functions R iz {q) fix normalization of the Baxter op- 
erator, while the residue functions E ± {q) define the en- 
ergy of the system (see Eq. ( JTc| ) below). 

Important properties of the solutions to (Q) can 
be deduced from the symmetry of the Hamilto- 
nian (^) under the cyclic permutations of particles, 
P^ , (zi, 2jv— 1, zn) = ^(z 2 , zn, z\) and their mirror 
permutations M^(zi, Zfi—i, %n) = ^(zn> •••) z 2 , Z\) 0. 
The eigenstate of the Hamiltonian diagonalizes the oper- 
ator F and the corresponding eigenvalue defines its quasi- 
momentum 



(q) = 2nk/N (8) 



with k integer in virtue of F N = 1. The operator M acts 
on the integrals of motion as qk — > (—l) k qk and maps the 
eigenstate of the Hamiltonian into another one with the 
same energy but different set of the quantum numbers, 
Mf q = with -q = (q 2 , -q 3 , (-l) n ?„), leading to 

E N (q) =E N (-q), Q q (-u,-u)=e* e »^Q„ q (u,u). (9) 



Inserting (j?]) into 
Then, substituting 



we find that E ± (q) = -E^(-q). 
) into (|]) and using hermiticity of 



the Hamiltonian (||), we obtain the energy as fjj 

E N (q) = Re[E+(-q)+E+(q)] . (10) 

3. Determining the eigenvalues of the Baxter operator, 
Q(u, u), we use the following ansatz || 



Q{u,u) 



d z z 



Q(z,z), (11) 



where Q(z,z) is a single- valued function on the 2-dimen- 
sional z— plane with z = z*. The Baxter equation (0) 
leads to the N-th order differential equation on Q(z, z) 



Iz^zdzfz 1 -* + z- j (zd z ) 

N 

= Y, lk 1k(zd z f- k Q(z,z) 

k=2 



Z J 



2(zd 2 



Q(z,z) 
(12) 



is fixed by a similar equa- 

1-3*, 9k = q* k - 

has three regular singu- 



The z— dependence of Q(z, z) 
tion in the z— sector with j = 

The differential equation (|l^ 
lar points: z = 0, z = 1 and z — 00. The solutions to 
( |l2| ) around z = 00 can be obtained from those around 
z = 0, using the symmetry of ( |l2"| ) under the transforma- 
tion z — > 1/z and — » (— l) fe ?fc. The solutions to ( |l2"| ) at 
the vicinity of z = and z = 1 have the form Q ~ z a and 
~ (1 — z) 6 , respectively. Solving the indicial equations 
one finds that the exponents a are N— times degenerate, 
(a — 1 + s) N = 0, while 6 are given by b\ — Nj — h — 1, 
62 = -^Vj +h — 2 and 6 m = m — 3 for m = 3, . . . , N. Since 
the exponents differ from each other by an integer, this 
may lead to the terms ~ In z and ~ ln(l — z). This is the 
case at z = 0, whereas at z = 1 the solutions to (U2) do 
not contain terms ~ ln(l — z) unless h = (1 + ny/)j2. 

Following 0, we define the fundamental set of solu- 
tions to Eq. (|l2|), Qm (z), (m = 1, AT), around z = as 



E 



Crn~l u k+l(z)(lnz) r 



(13) 



with the coefficients = T(m)/(T(k + l)T(m - k)) 

chosen for later convenience. The functions u m (z) are 
given by the power series u m (z) = f + X^feLi 2 ' u m,k with 
the coefficients u m ,k — Um,k(Q) satisfying the three-term 
recurrence relations with respect to k. Defining similar 
set of the fundamental solutions in the z— sector, we 
construct the general solution for Q(z, z) around z = as 



Q(z,z) 



x 

E 

m.m— 1 



(0) 
■in m 



(14) 



with being the mixing matrix. The fundamental so- 
lutions Qm (z) and Q { °/ (z) have a nontrivial monodromy 
around z = z = 0. For Q(z,z — z*) to be single-valued 
at z = 0, the monodromy should cancel in the r.h.s. of 
(|i"4|). To satisfy this requirement, Cnm should vanish for 
m > N —n, and have the following form for m < N+l — n 



N—n—m+l 



C (0) 



E 

k=0 



( — 2) k Otk+n+m-l 

kl(n- l)!(m- 1)! 



(15) 
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with ot\, ...,ajsf being arbitrary complex parameters. To 
fix the normalization of the function Q(z, z) we choose 
a N = l. Substituting @ and @ into @ we find 



0(^,2) 



In 



JV-l 



(zz) In 2 (zz) 



(iV-1)! (JV-2)! 



CCJV-1 +• 



(16) 

where ellipses denote the terms subleading at 2 = 0. 

As follows from (|ll|), the small— 2: asymptotics of the 
function Q(z, z) is in one-to-one correspondence with 
the singularities of Q(u,u). Substituting ( |l6|) into (|ll| ) 
and performing integration over |z| <C 1 one finds that, 
in agreement with ([7]), Q(u,u) has the iVth order pole 
i(j — 1) and uf = i(j — 1) 



at Ui 



Q(uf+e,uf+e) 



(ie) 



N 



{l+iea N _ x + 0(e 2 )}. (17) 



Matching this expression into (]?]), one obtains E + (q) = 
ajv_i(g). This leads to the following remarkably simple 
expression for the energy ( pp| ) 



E N (q) = Re [a N -i(-q) + a N -i(q)} , 



(18) 



where q = {%} and —q = {(— l) fc g^}. Thus, the en- 
ergy is related to the small— z asymptotics of the function 
Q(z, z), Eq. (0), or equivalently to the matrix elements 
of the mixing matrix ( p"5| ) in the fundamental basis (13). 

The fundamental set of solutions to Eq. (|l2|) around 
2 = 1 is defined similarly to (O) as 



Q${z)=z 1 S{l-z) b <"v m (z) 



(19) 



with m = 1,...,N. Here, the functions v n (z) are given 
by the power series, v n (z) = 1 + X)fc°=i(l ~ z ) k v n ,k with 
the expansion coefficients v n> k satisfying the N— term 
recurrence relations with respect to fc. To fix uniquely 
their solutions, we impose the normalization conditions 
v m {z) = l + 0((l-z) N - m+1 ) for m > 3. As was already 
mentioned, at h=(l+rih)/2 the solutions Qi\{z) become 
degenerate and one of them, Qi (z) for rih > 0, has 
to be redefined by including the additional ln(l — z)— term 

Q«(z) =^-i(i_ z )W-(n h +3)/a 

x [(1 - 2)"- ln(l - z) u 2 (z) + . (20) 

Using (|l^) and similar expressions in the z— sector, one 
constructs the general solution to JI4 ) for Im h ^ that 
possesses a trivial monodromy around z = 1 as 

Q(z,z)=(3 1 q[ 1 \z)q[ 1 \z)+P 2 Q 2 1 \z)Q ( 2 \z) 



N 



+ E €W7 mffl Q,! , (^>Q (1) 'C (1) 'g (1) .(2i) 

m,m=3 

At /i = (l+n/i)/2, or equivalently Im h = 0, the first term 
in the r.h.s. of (Ell) looks differently in virtue of (pOl) 



Q{z,2) = fa Q^(z)Q 2 1 \z) + Q 2 1) (z)Q^(z) +..., (22) 



where ellipses denote the remaining terms. 



Eqs. ([l4]) and ( pl| ) provide the solutions to the Baxter 
equation at the vicinity of z — and 2=1, respec- 
tively. Choosing z to be inside the region of convergence 
of the both expressions, we define the transition matrices 



N 



Q^(z) = J2Q nm Qg\z) 



(23) 



and similar for f2. For the fundamental set of solu- 
tions, Eqs. (13) and ([l9|), these matrices are uniquely 
fixed. They can be calculated as Q = [W^] -1 , with 



d n Qk ) i z o) and z being some reference point 



The transition matrices allow us to analytically con- 
tinue the solutions Jl4| ) from z = to z = 1. Substitut- 
ing ( p3| ) into ( p^) and matching the resulting expression 
for the function Q(z,z) into (Ell), we find the relation 
between the mixing matrices C^* and 



c (1) = [n{ q )f c m n($) 



(24) 



This matrix equation provides the quantization condi- 
tions for the integrals of motion g and fixes uniquely the 
eigenvalues of the Baxter operator. Replacing the mix- 
ing matrices by their expressions, Eqs. (li}) and (pl|), we 
obtain from (|24j) the system of N 2 equations involving 
(N - 1) parameters a n inside C (0) , 2+ (iV-2) 2 parame- 
ters /3i,2 and 7 mlft inside C*W, as well as (N — 2) integrals 
of motion qs,...,q^ (we recall that ffc = q^). Thus, the 
system ( pi| ) is overcomplete. Not only it allows to de- 
termine all parameters including the quantized q, but it 
also provides (2N — 3) nontrivial consistency conditions. 

4. Following this approach, it becomes straightfor- 
ward to solve the Schrodinger equation (^) for arbitrary 
SL(2, C) spins (j, j) and any finite N. In what follows, we 
shall present the results of our calculations of the ground 
state energy of (|l]) at j = and j = 1. Its value de- 
fines the energy of the N— reggeon states, En, entering 
Eq. (0). The detailed description of the full spectrum 
will be presented in the forthcoming publication. 

At N = 2, the equation ( |l2| ) can be solved exactly and 
the general expression for the function Q(z, 2) be writ- 
ten in terms of the Legendre Q— functions of the 2nd kind 



Q(z,z) = 



(1-z} 



■[Qh{z)Q- h {z)-Q X - h {z)Q 1 _- h [z)], 



(25) 

with Qh(z) = Q-h ((1 + 2)/(l — z)). Its expansion at 
z = and 2 = 1 matches Eqs. (|l6|) and (^l|), respectively, 
and leads to a\{q) = a\{— q) = E 2 (h)/2 with 

E 2 (h) = 4 Re h/>(/i) + </>(! - h) - 2ip(l)] (26) 



and e l6 * 2 = (— l) nh , in agreement with the known exact 
expression M. The ground state is unique and it corre- 
sponds to h= 1/2, E 2 = -16 In 2 and 6 2 = 0. 

For TV > 3, the analytical solution to (|lj) is not avail- 
able and one has to rely on the power series solutions, 
Eqs. ( |l3| ) and ([l9]). Performing thorough analysis of the 
quantization conditions (p4j) at N = 3, 4, 5, 6 with a help 
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<72 


iqz 


94 




96 


-E N /4 


o-]v/4 


N 


= 2 


.25 










2.77259 


16.8288 


N 


= 3 


.25 


.20526 








-.24717 


.9082 


N 


= 4 


.25 





.15359 






.67416 


1.3176 


N 


= 5 


.25 


.26768 


.03945 


.06024 




-.12751 


.4928 


N 


= 6 


.25 





.28182 





.07049 


.39458 


.5634 



TABLE I: Quantum numbers, q, and the energy, En, of the 
N— reggeon states in multi-colour QCD. The parameter crjv 
defines the energy of the lowest excited states, Eq. (p7|). 



of the method described in 1 10 (see Appendix C), we ob- 
tained the results summarized in the Table |. At N = 3 
they are in agreement with Q. 

The quantized values of the integrals of motion have a 
rich structure - they form the family of one-dimensional 
continuous trajectories in the (N— l)-dimensional space of 
q=(q2, <1n)- The "proper time" along each trajectory 
is defined by the imaginary part, Vh, of the conformal 
spin h = (1 + Ufc)/2 + ii>h- The total number of trajec- 
tories is infinite. Each of them is labelled by the set of 
integers including the quasimomentum 9n = 2irk/N and 
the parameter that defines the total 2-dimensional 
angular momentum of the N— reggeon state. In an 
agreement with our physical intuition, the ground state 
of the system belongs to the trajectory with On = and 
rih = 0. The energy along this trajectory is a continuous 
function of v% and it approaches its minimal value 
(=ground state energy En) at = 0. At the vicinity of 
Vh = 0, one finds an accumulation of the energy levels, 

E N (y h ) - E N = a N \v h \ d " , (27) 

with (In — 2 and the coefficients un given in the Table @. 

The ground state of the A^— reggeon system has dif- 
ferent properties for even and odd N . For even N , the 
ground state is unique and it belongs to the trajectory, 
on which the quantized q^,...,qN satisfy the relations 
Im g2fe+2 = q 2k+ i = with k = 1, (N - 2)/2. For 
odd N, the ground state is double degenerate. The cor- 
responding eigenstates belong to two different trajecto- 
ries, on which the quantum numbers satisfy the rela- 
tions Im<j 2 fc+2 = Reg 2 /c+i = and Img 2 /c+i ^ with 
k = 1, (N— 1)/2. In distinction with even N, quan- 
tized 53, 95, ... do not vanish and take pure imaginary val- 
ues on the ground state trajectories. Two ground state 
trajectories differ from each other only by the sign of 
<?2fc+i and the degeneracy of the ground state occurs due 
to the symmetry of the energy under q 2 k+i — * — <Z2fc+i, 
Eq. (||). The reason for this is M that at odd A" there 
exist two mutually orthogonal ground states, ^/ q ±} (zi) = 



q {zi)±^ - q (zi)]/2, which are invariant under the cyclic 
permutations, 8n — 0, and possess a definite parity un- 
der the mirror permutations, M^ q ±} (zi) = ±' ( £ q ±) (zi). In 
contrast, at even N one has ^- q (zi) = ^ q (zi) due to 
92/c+i =0 and, as a consequence, \£g — ^> q and 1 i''f > = 0. 

In high-energy QCD, in virtue of the Bose symmetry, 
the compound states of A" reggeized gluons has to be 
symmetric under interchange of any pair of reggeon co- 
ordinates and their colour indices. In the multi-colour 
limit, this symmetry is reduced to the invariance of the 
wave function under the cyclic and mirror permutations. 
The fact that the wave function of the state is factorized, 
as N c -~ xx), into the product of the colour tensor and the 
scalar function ^^^(^i, zn) implies that the both fac- 
tors have to possess the same parity under the cyclic and 
mirror permutations. This allows to distinguish the de- 
generate ground states according to their C— parity. The 
N— gluon states possessing the parity M = (— \) N have 
the same C— parity as the BFKL Pomeron, C = 1. For 
odd N, the ground states with the parity M = 1 have 
the C— parity of the Odderon state ||, C = — 1. 

Solving the system (p4|), we find the mixing coefficients 
and calculate the ground state energy from (p8J) . The re- 
sults are presented in the last two columns in Table |[ For 
En < (En > 0) the contribution of the N— reggeon 
state to the scattering amplitude A(s,t)/s, Eq.(|2|), in- 
creases (decreases) at high energy s. Our results indicate 
that for even N the ground state energy, En, is negative 
and it increases with A^ approaching the value -E200 = 
from below. For odd N, the ground state energy is posi- 
tive and it decreases with A^ approaching the same value 
E20C+1 = from above. This suggests that the ground 
state has the properties of an antiferromagnet. 

We conclude that, in the multi-colour limit, in the 
Pomeron sector, only states with even Af provide the 
contribution to (^|) rising with the energy s. Their in- 
tercept, however, is much smaller than that of the BFKL 
Pomeron. In the Odderon sector, the corrections to the 
"bare" Odderon A^ = 3 state due to higher states built 
from odd number of reggeons, N, become important. 
They increase the effective value of the Odderon inter- 
cept and lead to a a = 1 . Finally, in the both sectors, the 
contribution of the N— reggeon states to the scattering 
amplitude ceases to depend on the energy s as N — ► 00. 
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